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All references are textbooks and published papers 
Reference(s) used listed on each slide, lower left, as ref# 


NASA Themes and Libration Orbits 

Enterprises of Space Sciences (SSE) and Earth Sciences 
l\ are a combination of several programs and themes 
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• The Triana mission is the lone ESE mission not orbiting the Earth. 

• A major challenge is formation flying components of Constellation-X, MAXIM, TPE, and 





Earth Science Launches 
Low Earth Orbit Formations 
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SPECS (Submillimeter Probe of the Evolution of Cosmic Structure): 
Interferometer 1 km at L2 


Future Mission Challenges 
Considering science and operations 
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> Solar sail applications 

> Continuous control to reference trajectories 

> Tethered missions 

> Human exploration 
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# 3-7 


Gaussian Lagrange Planetary Equations 
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Potential Accelerations 
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Contains lag heating terms 
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#3-7 


Solar Radiation Pressure Acceleration 
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Other Perturbations 
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Inertial pointing spacecraft could have drastic changes in B c over an orbit 


Ballistic Effects 

Varying the mass to area yields different decay rates 

Sample: 100kg with area of 1, 10, 25, and 50m 2 , C d =2.2 



satlOOtoI.EIapsedDays (Days) 
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Matlab ODE 4/5 is a variable step RK 






ECEF - Navigation 
Topographic - Ground station 





Sr(t) = f(r,v)-f(r , 


Describing Motion Near a Known Orbit 
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This is important since it will be our starting point for everything 
that follows 



- Hill’s - RTN 

- Clohessy-Wiltshire - RAC 

- LYLH - RIC 


Describing Motion Near a Known Orbit 
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# 5,6,7 


Transforming the EOM 





Transforming the EOM yields 
Clohessy-Wiltshire Equations 
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a . (associated with T) has twice 

cos {nt) + — sin (nt) the amplitude of the x motion 



Relative Motion 

A numerical simulation using RK8/9 and point mass 
Effect of Velocity (1 m/s) or Position(l m) Difference 



- 0.045 - 0.040 - 0.035 - 0.030 - 0.025 - 0.020 - 0.015 - 0,010 - 0.005 0.000 

main.AlongArcTrackSeparatlon (Km) 
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Graphics Ref: Collins, Meissinger, and Bell, Small Orbit Transfer Vehicle (OTV) for On-Orbit 
# 7,9 Satellite Servicing and Resupply, 1 5th USU Small Satellite Conference, 200 1 






Consider two trajectories r(t) and R(t). 



# 5 , 6,7 


Lambert Problem 
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‘single pass’. 









EO-1 GSFC Formation Flying 

New Millennium Requirements 



re re 



The onboard flight software, called the EFF, shall provide the interface between the ACS / 
C&DH and the AutoCon™ system for Autonomy for transfer of all data and tables. 




EO-1 Formation Flying Algorithm 
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where F and G are found from Gauss problem and the f _& 9 
series and C found through universal variable formulation 




EO-1 AutoCon™ Functional Description 



- Performance Objectives 
• Flight Constraints 















Burn Decision and Planning 
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Inclination Maneuver Validation : Computed AV at node crossing, of 
~ 24 cm/s (114 sec duration), Ground validation gave same results 



Difference in EO-1 Onboard 
and Ground Maneuver Three-Axis AVs 



r# 10,11 






Formation Data from Definitive Navigation 

Solutions 
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separation over all 
formation maneuvers 
maintained to 3 km 


Formation Data from Definitive Navigation 

Solutions 
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Formation Data from Definitive Navigation 
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Enables Autonomous StationKeeping, 
Formation Flying and Multiple Spacecraft Missions 







NASA Libration Missions 
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• MAXIM Lissajous, Formation Flying of Multiple S/C 

• TPF Lissajous, Formation Flying of Multiple S/C 

(Previous missions marked in blue) 







# 1 , 12,13 Farquhar et al [1985] Trajectories and Orbital Maneuvers for the 1SEE-3/1CE, Comet Mission, JAS 33, No. 





Goddard Space Flight Center 



350,000km, Az~ 250,000km, Class I 
Spacecraft: Mass=1254kg, Spin Stabilized, 

Notable: First Ever Multiple Gravity Assist Towards LI 



Goddard Space Flight Center 



Spacecraft: Mass=818kg, Three Axis Stabilized, 

Notable: First Gravity Assisted Constrained L2 Lissajous Orbit; Map-earth Vector 

Remains Between 0.5° and 10° off the Sun-earth Vector to Satisfy 
Communications Requirements While Avoiding Eclipses 52 
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Aspherical Configurations (Position & Rates) 
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CR3BP Analysis of Phase Space Eigenstructure 

Near Halo Orbit 





Natural Formations: 

Quasi-Periodic Relative Orbits — > 2-D Torus 






Floquet Controller 

(Remove Unstable + 2 of the 4 Center Modes) 
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Deployment into Torus 
(Remove Modes 1, 5, and 6) 
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Deployment into Natural Orbits 
(Remove Modes 1, 3, and 4) 




O 

dJ 

GO 

6 S 




# 15 



Natural Formations: 
Nearly Periodic Relative Motion 
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Natural Formations: 
Expanding Vertical Orbits 



Chief S/C 





Fixed Relative Distance and Orientation 
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Others 




Formations Fixed in the 
Rotating Frame 



y (10 km) 



2-S/C Formation Model 
le Sun-Earth-Moon System 



Ephemeris System = Sun+Earth+Moon Ephemeris 










Linear Targeter 







Achievable Accuracy 
via Targeter Scheme 



Formation Distance (meters) 


Continuous Control: 

LQR vs. Input Feedback Linearization 
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— » Substitute x° (/) into EOMs and solve for u d ° (t) 
-+u d {t) = u;(t)-K(t){x{t)-T(t)) 

-> Apply the computed control input to the perturbed EOMs 
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Output Feedback Linearization 
(Radial Distance Control) 



Scalar Nonlinear Constraint on Control Inputs 

h(7(t),¥(t)^-u(t) T r(t) = 0 



Output Feedback Linearization (OFL) 
(Radial Distance Control in the Eohemeris Model ') 


00 



Critically damped output response achieved in all cases 
Total AV can vary significantly for these four controllers 



OFL Control of Spherical Formations 
in the Ephemeris Model 



Relative Dynamics as Observed in the Inertial Frame 





OFL Controlled Response of Deputy S/C 
Radial Distance + Rotation Rate Tracking 
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OFL Controlled Response of Deputy S/C 
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OFL Control of Spherical Formations 
Radial Dist. + Rotation Rate 
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Radial Distance (km) 



Inertially Fixed Formations 
in the Ephemeris Model 



Nominal Formation Keeping Cost 
(Configurations Fixed in the RLP Frame) 




Max./Min. Cost Formations 
(Configurations Fixed in the RLP Frame) 



Nomina! Relative Dynamics in the Synodic Rotating Frame 




Formation Keeping Cost Variation 
Along the SEM and L 2 Halo Families 
(Configurations Fixed in the RLP Frame) 
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x (10 km) A (to 6 km) 


Continuous Control in the Ephemeris Model: 
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Nearly periodic & quasi-periodic formations in the RLP frame 
Floquet controller: numerically ID solutions + stable manifolds 


Some Examples from Simulations 
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Periodic Reference Orbit 


h' 


+ 



+ 


&V 

cn 


O 


O 


$ 

• tH 

C/3 

§ 

& 

X 

§ 

1 


II 

II 

-H 



CD 

<n , *— i 

B § § 

• rH O 

•a, a a 


rx 3 ^ 

& a" 52 

S oj 2 
< * £ 


II II II 

^ 3 ^ 




.s ° 

S O 


I ^ 

II II 

PS N 









Centralized LQR Design 
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Centralized LQR Design 
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Disturbance Accommodation Model 
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Disturbance Accommodation Model 




r- 




#16,17, 20 




Motion of Formation Flyer With Respect to 
Reference Spacecraft, in Local (S/C-l) Coordinates 
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Stellar Imager Concept 

(Using conceptual distances and control requirements 
to analyze formation possibilities) 



ts M 


S £ 


> OQ 

<D 

c\ 




spacecraft, and environmental conditions 

#16,17, 20 




CD 

cn 

ru 

E 


_fD 

ID 

4-J 

CO 


£ 

QJ 

H 

4D 

O 

*h 


60 

• S QJ 

?H 1—1 

£ 2 

60^ 


O 


o 

4-> 

£ 

O 

u 

£ 


?H 

T5 

£ 

as 

v 

£ 

O 


43 

!h 

O 

£ 

as 

C/3 

£ 

•i-H 

as 


£ 

O 

T2 OJ 
as e£ 

3 H 

r-* 

£ 

C/3 2 

•i-H 

aj 

x JS 

£ 


•i—l 


° £ 

bo *53 

£ QJ 


o 
• 1— < 
-4-* 

•4— • 

as 

Jh 
•1— ( 

as 

S 

•rH 

44 

4-> 

» 

as 

as 

g 

s 

Vh 

•T— ( 
60 

M-l 

O 

£ 

o 

O 

Mh 

£ 

QJ 

46 

T5 

£ 

CJ 

Ph 

Mh 

QJ 

X 

QJ 

QJ 


i— i 

cn 

44 

4— > 

QJ 

44 

QJ 



QJ 

q. u 

^ QJ 

-M 

8 3 

U g 

> « 
< S 

£ bo 

-2 .S 

_aj 

03 





Stellar Imager Mission Study Example Requirements 
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Results - Libration Orbit Maintenance 
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Represents both 0.5 and 4 km focal lengths 












Results - Formation Slewing 



Hub estimation 0.5 km Drone estimation 0.5 km 

separation / 90 deg slew separation / 90 deg slew 
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Results - Formation Reorientation 
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Results - Formation Reorientation 
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Summary 

(using example requirements and constraints) 
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tumroi; ior ai couia De accomplished with better sensors to 
lessen the effect of the process and measurement noise. 


Summary 

(using example requirements and constraints) 
the controller and varying the maneuver intervals 
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• system reliability and lault detection 

• Nonlinear control and estimation 

• Second and third control tiers and new control 
strategies and algorithms. 
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DRO Mission Metrics 
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Distant Retrograde Orbit (DRO) 



Asymptotic Orbit in 
the Restricted Three- 
Body Problem 


Earth Distant Retrograde Orbit (DRO) Orbit 
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DRO Formation Control Analysis 
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DRO Formation Control Analysis 
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DRO Formation Control Analysis 
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Formation Control Analysis 



03 



d) Resize from 25km to 50km 

e) Maintain 50km sphere (strict PD control) one month 

f) Maintain 50km sphere (loose control) one month 
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DRO Formation Control Analysis 

> Earth/Moon L4 Libration Orbit 

> Spacecraft controlled to maintain only relative separations 
^ Plots show formation position and drift (sphere represent 2 
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General Theory of Decentralized Control 



•Data transmission requirements are ‘ W 

MINIMIZED. Node 4 
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DST/Numerical Comparisons 
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Libration Point Trajectory Generation Process 










irajectory ueneratec 
with Manifold States 




JWST DST Perspective 





1 raj ectory Lrenerated starting \ 

with Manifold States Starting Point 
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Coordinate system transformation 
Orbital to inertial coordinates 






Principles Behind Decentralized Control 







The LQG Decentralized Controller Overview 



TO NETWORK 




